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ABSTRACT 

We present a generic algorithm for performing astronomical image registra- 
tion and pointing refinement. The method is based on matching the positions 
and fluxes of available point sources in image overlap regions. This information 
is used to compute a set of image offset corrections by globally minimizing a 
weighted sum of all matched point-source positional differences in a pre-specified 
reference image frame. A fast linear sparse matrix solver is used for the mini- 
mization. From these corrections, the pointings and orientations of images can 
be refined in either a relative sense where pointings become fixed (registered) 
relative to a single input image, or, in an absolute sense (in the ICRS) if absolute 
point source information is known. The latter provides absolute pointing refine- 
ment to an accuracy depending on the robustness of point source extractions, 
match statistics, and accuracy of the astrometric catalog used. The software is 
currently used in the Spitzer image processing pipelines, although it is adaptable 
to any astronomical imaging system which uses the FITS image format and WCS 
pointing standard. 

We test the algorithm using Monte Carlo simulations and compare them to 
image data acquired with the Infrared Array Camera (IRAC) on board the Spitzer 
Space Telescope. We find that dispersions in matched source separations after 
refinement are entirely consistent with centroiding errors in source extractions, 
implying that systematic uncertainties due to inaccurately calibrated distortions 
are negligible. For these data, we predict refinements to better than ~ 70 and 
~ 280 mas {2a radial) for the IRAC 3.6 and 8/xm bands respectively. These 
bands bracket two extremes in available source matches and, for the data under 
study, correspond to an average of about 55 and 8 matches per image in these 
two bands respectively. 
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Subject headings: methods: data analysis — techniques: image processing — catalogs— 
astrometry 

1. Introduction 

The need to register or co-align images in astronomy to an accuracy better than the 
Nyquist sampling density of a detector's point response function is pivotal if one wants 
to maximize the resolution and signal-to-noise attainable in a mosaic of co-added images. 
The absolute accuracy is inhibited by instabilities in telescope pointing control, tracking 
sensors, and how these behave in the science instrument frame in the presence of thermo- 
mechanical disturbances. For instance, imaging detectors on both the Hubble 1 and Spitzer 2 
Space Telescopes provide an absolute pointing accuracy of ~ 0.5-1 arcsec (1-cr radial). For 
comparison, the highest attainable resolutions 3 are ~ 0.02 and 1 arcsec for these telescopes 
resepectively. Without further refinement, the current pointing accuracies are insufficient to 
exploit the near-diffraction limited resolution capabilities the detectors can provide. Factors 
of 10-20 improvement in pointing are required for optimal image registration. 

Good image registration enables extraction and position determination of sources to 
fainter flux levels for a given signal-to-noise ratio. Comparison or registration with astro- 
metric sources whose positions are known to better than a few percent of the observed image 
pointing uncertainty also allows refinement of image frame pointings in the International 
Coordinate Reference System (ICRS). Absolute pointing refinement can be achieved to an 
accuracy approaching that of the astrometric catalog used or better given good match statis- 
tics. This alleviates possible ambiguities when performing cross-identification/correlation of 
extracted sources across wavelength dependent catalogs. Furthermore, the accurate place- 
ment of slits for follow-up spectroscopic studies requires source positions known to better 
than a few tenths of an arcsecond in the ICRS, or until the desired positional accuracy is 
limited by centroiding error in the array frame. 

Broadly speaking, image registration methods can be loosely divided into three classes: 
algorithms which use information in pixel space directly, e.g., by correlating common objects 
(Barnea and Silverman 1972); algorithms which attempt to match features or identified parts 
of objects known as "graph-theoretic" methods (Brown 1992); and algorithms which use the 

1 http://www.stsci.edu/hst/observatory/pointing 
2 http://ssc. spitzer. caltech.edu/documents/SOM/ 

3 Measured in terms of the Full- Width Half-Maximum (FWHM) of the center of an Airy disk pattern 
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frequency domain e.g., methods based on computing cross-correlation power spectra via the 
Fast Fourier Transform (Kuglin and Hines 1975). The two conventional methods for regis- 
tering images in astronomy involve either interactively identifying common point sources in 
overlapping image fields (i.e., object correlation), or, using actual detector acquired point- 
ings (with inherent uncertainties) directly and estimating relative image offsets therefrom. 
These methods are available in most data reduction packages (e.g. IRAF - Davis (1996); 
STARLINK - Bly et al. (2002)) and are mostly limited by accuracies in point source cen- 
troids, sufficient match statistics, telescope pointing control, or, subject to random-walk 
(cumulative) uncertainties. 

A number of robust methods for astrometric calibration of single images have also been 
implemented in the commonly used data reduction packages (Veran and Rigaut 1998; Valdes 
1998; Bustos and Calderon 2000), although an automated, self-consistent means for simul- 
taneous registration and refinement of multiple astronomical images comprising a mosaicked 
region is generally lacking. The single image methods assume one has sufficient numbers of 
astrometric matches to mitigate against uncertainties in the centroids of extracted sources. 
We have gone a step further by using the available point-source content to obtain a global re- 
fined solution for all frames such that all matched point-source positional discrepancies in all 
frame-to-frame and frame-to-absolute overlaps are minimized. Combining the relative and 
absolute source information reduces the demand on having sufficient astrometric matches. 
This becomes important for mid to far-infrared imagery where cross-wavelength astrometric 
calibration is often unreliable due to differing sensitivities, source populations and detector 
point spread function sizes. 

We have developed an algorithm to register and refine simultaneously the pointing 
of an ensemble of astronomical images to accuracies better than that inherent in point 
source centroid uncertainties (and dictated by point source match statistics). This paper 
describes the global minimization algorithm and presents a case study using data from 
Spitzer's Infrared Array Camera (IRAC). The outline is as follows. Section 2 describes the 
algorithm and pointing refinement accuracies expected therefrom. Section 3 validates the 
algorithm using a Monte Carlo simulation of IRAC data. Section 4 compares these results 
to real observations acquired with IRAC. Discussion and conclusions are given in Section 5. 



2. Algorithm 

The algorithm has been implemented into a stand-alone software package called point- 
ingrefine. The main inputs to pointingrefine are a FITS image list, with each FITS image 
containing the standard World Coordinate System (WCS) keywords (Greisen and Calabretta 
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2002; Calabretta and Greisen 2002), an accompanying list of flux-calibrated point-source 
extraction tables, optional astrometric source list, and configuration parameters. The soft- 
ware expects point source extraction tables adhering to the format generated by the Spitzer 
Science Center (SSC) source extractor. The complete pointingrefine package is available 
at http://ssc.spitzer.caltech.edu/postbcd/. This includes programs to perform point source 
extraction as well, although these will not be described here. A general overview of the 
processing steps involved in the global minimization algorithm (the heart of pointingrefine) 
is shown in Figure 1. In the following sections, we expand on some of the more important 
computational steps of this algorithm. 



2.1. Set-up and Point Source Matching 

Prior to source matching, we first compute the total number of image-pairs expected 
to be overlapping (which could potentially contain correlated/common sources) in the input 
ensemble of images. This is used for efficient a-priori memory allocation. Given a number 
of images, N imgs , the maximum number of distict frame pairs that can mutually overlap is 

■N maxpairs ^^irngs i^-^irngs 1) • (1) 

This maximum occurs when all images are stacked more-or-less on top of each other. For a 
sparse mosaic, this number is smaller and thus puts less of a burden on the required memory. 
The total number of potential overlapping frame pairs is found by computing the distances 
between the centers of every image pair and finding whether the distance is less than the 
typical radius of a circle enscribing an image. 

Source positions from all extraction tables (including absolute astrometric references 
if specified) are sorted in declination. This preconditioning speeds up the source matching 
procedure, converting it from an 0(N 2 ) to an approximate O(N) linear process. Every pos- 
sible overlapping image pair is scanned for common point sources in the RA, Dec coordinate 
system. Both position and flux matching is performed. The position match step attempts 
to find sources falling within a nominal search radius (typically several times the root sum 
squares of prior image pointing uncertainty and extraction centroid error). If more than one 
match is found within the search radius, no match is declared, due to possible ambiguity. In 
other words, only singly matched sources within the search radius are used. Sources are also 
simultaneously matched in flux. A flux match is satisfied if two fluxes fall within a maximum 
tolerable relative flux difference threshold. The software includes options to re-scale the in- 
put fluxes of astrometric references to be commensurate with those of actual extractions. A 
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minimum of two matches per image is enforced to declare a correlated image pair, since this 
is the minimum required to estimate a rotational offset between the pair unambiguously. 

In preparation for the global minimization step (see below), all point source match 
positions and uncertainties are mapped into the pixel coordinates of a Cartesian reference 
image frame. The definition of this reference depends on whether absolute or relative re- 
finement is desired. The pointingrefine software distinguishes between these options from its 
given inputs. Two inputs are required for absolute refinement. First, the software expects 
a Fiducial Image Frame (FIF). This is a file listing various WCS parameters that define the 
tangent point and boundary of a "fiducial image" encompassing all images. Second, a list 
of astrometric (absolute) sources that fall within the FIF boundary is required. The FIF 
acts as an effective input image whose sources are the astrometric references. Refinement 
with respect to such a FIF ensures both absolute and relative refinement amongst images. 
In relative refinement mode, all input images can be registered and refined with respect to 
a single input image. In this mode, either the software selects an image from the input list 
which is maximally correlated (has most overlaps) with other images, or the user can specify 
their own. 



2.2. Global Minimization 

Consider the simple three-image mosaic in Figure 2. Image "1" defines the "fiducial" 
reference frame. The circles represent point sources detected from each overlapping image 
pair transformed into the reference frame of image 1. The filled circles are sources extracted 
from image n and the open circles are sources extracted from either image 1 or m. The 
matched source pairs are shown offset from each other to mimic the presence of random 
pointing uncertainty in each input image. These are the offsets we wish to compute and use 
to correct each frame pointing. A one-dimensional representation of the projection geometry 
in the reference image frame is depicted in Figure 3. All input images have sizes in the 
reference frame that depend on their distance from the reference image tangent point. The 
projected linear size scales with angular distance 6 t as ~ 1 + tan 2 ^. The pointingrefine 
algorithm appropriately accounts for the inflation of centroid uncertainties and separations 
of correlated sources when projected into the reference frame. A potential problem for the 
algorithm is if image sizes and mosaic extents are large enough to cause a non-uniform 
dependence in scale over the region in which an input image is projected. The dependence 
of this effect with 9 t , however, is weak. For instance, the projected scale varies by < 1% 
across an image 30' in size at angular distances 9 t < 30°. 

The main assumption of the algorithm is that the random uncertainty in measured twist 
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angle 4 of an individual image frame (59 m in Figure 2) is small enough to ensure sin 59 m w 59 m 
(see below). 50 < 1' is a good working measure for the purposes of this algorithm. This is 
justified for the Spitzer science instruments where the absolute twist angle uncertainty 5 is 
typically < 30" (1-<t). 

In a rectilinear coordinate system (say that defined by the frame of image 1 in Fig- 
ure 2), the positions of a point source i detected from images m and n are related via the 
transformation: 



X'- 

y\ 



x\ 

y r i 



X 



+ 



cos 59 m - sin 59 m 
sin59 m cos 59 m 



+ 



5X m 
5Y rn 



(2) 



where 59 m is a rotational offset, 5X m , 5Y m orthogonal translations, and (x™ , y™) the center 
coordinates of image m in the reference image frame. The rotation 59 m is measured in a 
counterclockwise sense so that a rotation followed by orthogonal translations in x and y of 
image m will align the sources (open circles) detected therein with those detected in image 
n. 

The assumption of small twist angle uncertainty (see above) implies sin 59 m « 59 rn and 
cos 59 m ~ 1 and thus, the pair of equations defined by (2) can be linearized in 59 m as follows. 



x ™^~ x n = xf _ {yT _ yT)59 m + 5x m 

(3) 

yT-^y? = y? + (xT-x™)59 m + 5Y m , 

where (x™, y") represents a position, corrected for unknown image rotation and orthogonal 
translations. 

We define a cost function L, representing an inverse-variance weighted sum of squares 
of all matched point-source positional differences from all overlapping image pairs (m, n): 



4 We define the "twist angle" as the relative orientation of an image in an orthogonal coordinate system, 
not the conventional position angle measured in the ICRS 

5 From Jet Propulsion Laboratory Internal Document: SIRTF Instrument Pointing Frame Kalman Filter 
Execution Summary (IPF Team), report ID01M095, October 2003 
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L ~ ^] \ a m,n [ X i ] + a m,n [Vi Hi ] f + ^apriori i (4) 

m,n(n>m) i L ^ > 

where 

Aaf'" = a\xT)+o 2 W) (5) 

and the a 2 represent variances in extracted point source centroids. The summations in 
Equation 4 are over all matched point-sources from all image overlaps, including sources 
matched to any astrometric references present in the FIF if available. The FIF or absolute 
reference frame (as defined in Section 2.1) is treated as an effective input image when source 
matching is performed. 

The additive term L apr i or i in Equation 4 represents an "a-priori" weighting function 
which makes use of actual measured pointing uncertainties of the images. This function is 
defined as: 



^ f(5X m ) 2 (5Y m ) 2 (56 m ) 2 (5X n ) 2 (5Y n ) 2 (56 n ) 2 \ 

Lapriori ~ ) ~~Z2 ' T2 1 3 ' 3 1 3 1 3 f ' W 

m ,n(n>m) L ° Xm ^ m C 6m ° Xn ° Yn °^ > 

where the u 2 m) a 2 n (j = X, Y, 8) represent measured pointing variances (in the ICRS) 
transformed into the reference image frame. The purpose of including L apriori is to avoid over- 
refining or biasing those images whose inherent measured pointing uncertainties are already 
small (within nominal requirements). In other words, those images whose uncertainties 
are known to be small a-priori, will have have a larger contribution to L apriori relative to 
the correlated source term (double sum in Equation 4). Consequently, the solution will be 
biased towards minimizing L apriori and not the correlated source term which could potentially 
degrade expected image offsets (and final refined pointings). In the limit L — > L apr i or i, the 
global minimum will be closer to 66 ~ 5X tv 5Y 0. Vice- versa, large pointing uncertainties 
will bias the solution towards the correlated source term where refinement via point-source 
matches is obviously needed. 

Equation 4 can be rewritten in terms of physical image offsets 58 m , 5X m and SY m for 
an arbitrary image m via use of Equations 3 and 6: 
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(7) 



The cost function denned by Equation 7 can be treated as a standard x 2 statistic to the 
extent that point source centroid uncertainties are independently random and Gaussian, a 
good approximation in this case. The probability density function for L about its global 
minimum is therefore the x 2 distribution for v degrees of freedom where 

v = 2N matches — S(N imgs — 1). (8) 

^matches is the total number of point-source matches in all image overlap regions and N imgs 
the total number of images containing the detected matches (including the reference image). 
The total number of parameters to solve is actually 3(N irngs — 1) since we have three offsets 
(56 m , 5X m and SY m ) for every image m and we exclude the reference image which by 
definition is constrained to have 69 = 0, SX = and 5Y = 0. 

Our aim is to minimize L with respect to all image offsets 59 m , 5X m and 5Y m for every 
correlated image m. At the global minimum of L, partial derivatives with respect to these 
three offsets for each image m are required to vanish: 



dL dL dL 

d~59^ ' d~5X^~ ' d~5Y^~ ( ) 

Evaluating these partial derivatives leads to a set of three simultaneous equations for each 
image in the ensemble. For N imgs , we therefore have 3(N imgs — 1) simultaneous equations 
in 3(N imgs — 1) unknowns. Each image m of an ensemble m = 1..M could be potentially 
correlated (contain sources in common) with any other image n where n = 1..N; m ^ n 
and M = N = N imgs — 1. From the conditions defined in Equation 9, the simultaneous 
system of equations used to solve for the offsets of all images m, can be represented as the 
matrix equation: 



- 9- 



r>m=l 
C m=l 



Am=l 



B 



Am=l 
Ay 



C m=l 



An=l 
Dn=l 
n=l 



at 1 

bt 1 



An=l 
Ay 



Cy =1 



\n=N 



r>n=N 



An=N 
A X 
r>n=N 



a 



n=N 



Ay =N \ 



(jn=N 



Am=M Ain=M Ain=M 

A e Ji X Ay 

6 X 

Qm=M q £jm=M 



( 59 m= 
5X mz 

5 ym-- 



1 \ 

=1 

= 1 



V 



^Qrn=M 



J 



1 \ 



,j,m=M 
^ A 

,T f m=M 
w B 

o 



/ 

(10) 



Equation 10 is of the form M.X = ^. Elements of the coefficient matrix M and the right- 
hand-side column vector ^ are given in Appendix A. The solution for the column vector 
X with 3(N imgs — 1) unknowns can be obtained using linear sparse matrix methods since, 
depending on the mosaic geometry, a large number of the matrix elements can be zero. In 
general, the fraction of zeros will be ^ § (> 22%). The ps 22% minimum will occur when 
every image of the ensemble is correlated with every other, such as in a coadded stack. The 
level of sparsity in M will increase with non-zero elements along a block diagonal if one 
desires to tie and refine images with respect to astrometric absolute references alone. In this 
specialized case, n = the reference image and all elements with superscript n in M will be 
zero. In general, the minimum and maximum possible number of zeros in the matrix M are 



N mm = 2(N mgs -l) 2 and (11) 

Nmax 9(-/Vj m g S 1) 7 \Nimgs 1) 

respectively. 

We use the UMFPACK 6 library to solve the matrix equation 10. This is adapted for 
solving large nonsymmetric matrix systems. The library includes an iterative scheme to cor- 
rect solutions for the inevitable accumulation in round-off error during the LU-factorization 
stage. 

We also compute the full error- covariance matrix for all image offsets (59 rn , 5X m , 5Y m ), 
which is given by the inverse of the coefficient matrix: C = M -1 (e.g., Press et al. (1999)). 
Variances in each offset are along the diagonal of C and covariances are given by off-diagonal 



6 htt p : / / www . cisc . ufl . edu /research/ sparse / umfpack / 
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elements. The covariance matrix C is computed using the same matrix solver as above by 
repeatedly solving for each unknown column X c in M _1 such that M.X C = I c where I c is 
the corresponding column in the identity matrix. An analysis of the covariance matrix using 
real data is presented in Appendix B. 

Images which are not correlated (lack point source matches) with any others in the 
input ensemble cannot contribute to the globally minimized cost function L. For these 
images, the pointingrefine software explicitly sets their reference frame offsets to zero and 
no refinement of their pointing is possible. As a further detail, the inclusion of the a- 
priori weighting function L apr i or i (Equation 6) guarantees that the matrix in Equation 10 
will be non-singular. The priors provide at least one measurement per image. If these are 
omitted from the cost function (Equation 4), then there are cases where the matrix could 
be singular. This can occur if the input image ensemble contains clusters of correlated 
images disjoint from each other with a non-contiguous/broken path existing between the 
clusters. This situation leads to an under-represention of images across the full simultaneous 
system of equations and the determinant will be zero. As indicated in the processing flow of 
Figure 1, if this occurs in absolute refinement mode, a second pass computation is attempted 
and only those frames which contain absolute astrometric matches are used. No relative 
frame-to-frame information is used and images are refined based on their absolute source 
content alone. The matrix effectively becomes block diagonal and a non-zero determinant is 
guaranteed. In relative refinement mode, no attempt is made to perform registration within 
each disjoint sub-ensemble. Instead, the software will abort with a message indicating that 
disjoint clusters exist. 

Once image offset corrections (59 rn , 5X m , 5Y m ) and associated uncertainties have been 
determined in the reference image frame, the final step involves refining the celestial pointing 
and orientation of each image m. This is performed by correcting the pointing centers (x™, 
y™) of each image in the reference image frame via Equation 3, i.e., 



and then transforming back to the sky to yield refined pointings in the ICRS. Image ori- 
entations are refined in a similar manner. For these, we correct and transform two fiducial 
points per image to determine uniquely the refined position angle. The main outputs of 
the pointingrefine software are additional WCS keywords written to FITS image headers 
representing refined pointings and orientations on the sky (see processing flow in Figure 1). 



x™(new) 
y™(new) 




y?(old) + SY 



m 



(12) 
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2.3. Optimization and Expectations 



The accuracy in pointing refinement or registration can be severely limited by possible 
systematics. For example, inaccurately calibrated image scale and/or distortions known a- 
priori to be position dependent will bias extraction centroids and have adverse affects on 
the separations of bonafide source matches and hence final globally minimized solutions. A 
method to test for possible contamination from systematics involves examining distributions 
of matched source separations after refinement and comparing these with that contributed 
by (random) centroid uncertainties. This will be performed using real data in Section 4. 

The presence of absolute-astrometric point sources are an important ingredient for re- 
finement and registration in general. These reduce the potential for a biased random walk 
in refined pointing with distance of an image from a fiducial reference. Such effects are also 
alleviated by using appropriate prior pointing uncertainty information in the global mini- 
mization cost function (see Equation 4). Absolute astrometric references provide a set of 
"anchor points" to which all extractions will be attracted. If one desires to perform absolute 
refinement with little or no astrometric references, then sufficient numbers of images are 
needed to increase the number of potential matches in frame overlap regions for a single 
point source. This is needed in order to approach a normal distribution about the expected 
absolute source position. In other words, large numbers of correlated source positions will 
ensure that the mean position of a correlated source cluster is close to the "truth" , or that 
which will result after refinement (assuming no systematics as discussed above). Due to the 
rarity of cases with numerous image overlaps providing good normally-distributed sampling, 
it is advised to use astrometric references wherever possible. 

The accuracy to which we can refine the pointings of an ensemble of mutually corre- 
lated images predominately depends on the number of point source matches available (both 
relative and absolute). With randomly distributed uncertainties in point source centroids, 
it is expected that the mean separation between matched source positions is approximately 
Gaussian after refinement, by virtue of the central limit theorem. In this limit, the (radial) 
uncertainty in refinement 7 of a single image will scale as 




(13) 



7 Approximated as the uncertainty in the mean source match separation with source positions weighted 
by their inverse variances. 
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where N ext and N a b s represent respectively the number of frame-to-frame and frame-to- 
absolute source match pairs in all overlap regions associated with the image, a ext is a typical 
source extraction centroid uncertainty, and a a ^ s is an astrometric source position uncertainty. 
This approximation assumes that a ext > or a a t> s > when either N ext > or N a i> s > 
respectively. For no astrometric matches, we set N a b s = and a a bs = 0, and Equation 13 
reduces to o ext \/ 2 j N ext . We expect to measure source extraction centroids to better than 
« 0.1 pixel, (k, 0.121" for Spitzer's IRAC focal plane arrays). If we assume for example 
astrometric positional errors of ~ 0.2" (conservatively speaking), then to refine image point- 
ings to an accuracy better than 0.1" will require at least five astrometric point-source 
matches per frame if N ext = 0, or less, if N ext > 0. One can see that with more point-source 
matches, the better the refinement. This assumes the observational setup allows for sufficient 
frame-to- frame overlap to ensure good numbers of relative matches N ext . If this is not the 
case, one will have to resort to using pure astrometric (absolute source) matches alone. 

To summarize, corrections for optimal refinement will effectively be given by the mag- 
nitude of frame pointing uncertainties with errors approximated by Equation 13. The latter 
assumes that in the limit of increasing number of matches, the mean source separation per 
image overlap region is approximately normally distributed with an uncertainty determined 
exclusively by point source centroid uncertainties. Any position-dependent systematic offset 
between source matches such as nonuniform pixel scale or inaccurately calibrated distortion 
will limit the refinement accuracy to the size of systematic error involved. 



3. Validation Against a Monte-Carlo Simulation 

We quantitatively assess the performance of the above algorithm using a simulation of 
1000 mosaicked images, with each image's coordinates modeled with an uncertainty drawn 
from a Gaussian distribution. The simulation is generic in that it represents a good overall 
representation of the type of data that could be acquired with modern optical/near-infrared 
detectors to moderately faint magnitudes {m opt ~ 23 or m nearIR ~ 19). To facilitate a 
comparison with real observations in Section 4, we have chosen to model the source count 
distribution and detector properties with that expected (and more or less observed) in the 
3.6/im band of Spitzer's IRAC instrument. A more detailed description will be given in 
Section 4. 

The "truth" source flux-density distribution was simulated using the models of Xu et 
al. (2003). These assume a high galactic latitude stellar model, several galaxy luminosity 
functions depending on galaxy morphological type, and exploit a large library of spectral 
energy distributions. This simulation was used extensively for predicting Spitzer source 
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populations (Lonsdale et al. 2003) and for validating processing pipelines. 
The simulation steps are as follows. 

1. A "truth" list of random source flux densities was generated using the models of Xu 
et al. (2003) covering an area of 0.45 square degrees down to a flux density of 10/x Jy 
at 3.6/im (equivalent to ~ 18.6 magnitudes in the Vega system). 

2. Truth sources were assigned both random and correlated positions within the area 
to be mosaicked. Galaxies were assumed to have a weak correlated component with 
amplitude (excess above random) based on an empirical if-band 2-point angular cor- 
relation function: w(6) ~ 0.001(#/deg)~ a8 . A pixel scale of 1.21" (characteristic of 
the IRAC arrays) was used when mapping sources to the pixel frame of the mosaic. 

3. Truth sources were convolved with a point spread function (PSF) 8 scaled by the appro- 
priate source flux. PSF-convolved truth sources were mapped into the mosaic frame 
with no pixel resampling. 

4. 1000 (256 x 256 pixel) "truth image" frames were generated from the full mosaic area 
using a dither and mapping strategy which assumed ~ 60% adjacent image overlap. 
The optimal map geometry for this 1000-image set was generated by D. Shupe (2002, 
SSC, private communication). RA, Dec, and Twist angle information is also attached 
to each image at this stage. 

5. A new set of 1000 images was generated (our simulated control sample) with uncertain- 
ties added to image pointings and position angles. Pointing offsets (prior uncertainties) 
were modeled as Gaussian random deviates along each independent orthogonal image 
axis. These were drawn from a zero-mean Gaussian distribution with a = 0.85" per 
axis. This choice for a is based on a pre-launch pointing knowledge of ~ 1.2" (1-cr ra- 
dial) for Spitzer in the science instrument frame 9 . Even though Spitzer can now actually 
do better than this by a factor of ~ 1.5 (Section 4), the ultimately refined pointings 
are independent of the magnitude of reasonable simulated prior uncertainty assumed 
(see below). Twist angle uncertainties were modeled as Gaussian with a = 20". A 
smoothly varying background adjusted with the expected Poisson and (IRAC band-1) 
read-noise per pixel was added to each image. 



8 Made from in-flight IRAC band-1 observations by the IRAC Instrument Support Team at the SSC. It 
has FWHM of - 1.66" and - 42% central pixel flux. 

9 Jet Propulsion Laboratory Internal Document: Operational Implications of the Time-dependent Pointing 
Behavior of SIRTF; C. R. Lawrence et al. Version 1.00; August 22, 2000 
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6. An absolute source list (representing astrometric references) was generated by taking 
the brightest truth sources which gave a mean density of ~ 50 sources per 5.2' x 5.2' 
image region. This resulted in 3030 truth sources. To imitate an astrometric catalog, 
the sources were assigned with positions modeled as truth ± u.a along each axis with 
u drawn from a Gaussian distribution with a = 0.06". This is typical for sources in 
the 2MASS point source catalog to K s ~ 15 (see Section 4 for details regarding the 
2MASS catalog). 

7. The SSC point source extractor was used on each simulated (control) image to extract 
sources above a threshold of ha. This resulted in m 40 extractions per frame. 

The pointingrefine software was executed on the 1000 image control sample. A source 
match radius of 3.5" was used to comfortably accomodate prior image pointing errors and 
extraction centroid errors (typically 0.15" 1-a per axis). Simultaneous flux matching was also 
applied between frame-to-frame and frame-to-absolute (astrometric) matches with maximum 
flux difference thresholds of 5 and 10% respectively. A zoomed-in (2.8' x 3.8') section of our 
1000-image simulation (with image pointing errors) is shown in the left panel of Figure 4. On 
the right is the same section after pointing refinement. The increase in resolution is dramatic. 
There is a factor of ~ 6 decrease in mean source match separation leading to more localized 
point source flux distributions and detectability to fainter levels. In this test, the surface 
brightness is increased by factors of ~ 2 — 3.5 for detected sources after refinement. This is 
as expected given that the 1-a radial image pointing uncertainty is of the order the input 
pixel size, and the instrinsic PSF has FWHM ~ 1.66". 

To get a more quantitative assessment of the performance of pointingrefine, we compare 
the distribution of separations between image center pointings of "truth" and simulated 
(control sample) images before and after refinement. This is shown in the top panel of 
Figure 5. Two different runs of pointingrefine were performed based on the number of 
(brightest) extractions used per image. One gave an average of ~ 10 relative (frame-to- 
frame) and ~ 20 absolute-astromeric source matches per image (the "10/20-match" case), 
and the second resulted in an average of ~ 2 relative and ~ 3 absolute matches per image 
(the "2/3- match" case). It should be noted that two matches per image is the absolute 
minimum to determine unambiguously two orthogonal shifts and an orientation per frame 
(i.e. the number of degrees of freedom (#dof) = 2N matc h es — 3 = 1). The dispersion in image 
separation relative to truth after refinement for the "10/20-match" case is ~ 65 mas 10 (1-a 
radial). For the "2/3-match" case, this is ~ 110 mas. For the given match statistics, these 



1 mas = 1 milli-arcsecond 
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numbers are more or less consistent with the simple scaling relation given by Equation 13. 
Offset distributions along each axis are shown in the lower panel of Figure 5 where the open 
circle represents the 2-a contour for the "2/3- match" case. 

By comparing differences in image radial separations before and after refinement, we 
find that over 890 images have their pointing refined to better than 95% relative to "truth" 
(for the "10/20-match" case). The top panel in Figure 6 shows the distribution in frac- 
tional refinement. This quantity is defined as the ratio of separations: 1 — D(refined — 
truth) I'D {unrefined — refined). For the "2/3- match" case, slightly less than half have the 
same amount of refinement, although most images are refined to better than 80%, corre- 
sponding to a discrepancy of ~ 180 mas within truth image positions. Thus the refinement 
is very good, even with minimal matches. 

The bottom panel in Figure 6 shows distributions in matched source radial separations 
before and after refinement. Uncertainties in both image pointing and source extraction 
centroids (or radial separation) contribute to the dispersion in the unrefined distribution. 
Image pointing uncertainties dominate, with a 1.2" contribution compared to ~ 0.2" for 
extraction centroids (both 1-er radial). After pointing refinement, the dispersion in radial 
source separation is expected to be dominated exclusively by extraction centroid errors and 
indeed, the distributions confirm this for both the "2/3" and "10/20-match" case (Figure 6 
bottom left). The distribution for the "2/3- match" case is narrower since such matches 
were performed using the brightest extractions per image, and these inherently have better 
determined centroids. 

The separation between two source extraction centroids (xi, i/i) and (xj, yj) is given by 

r = [{x l -x 3 f + {y l -y 3 f] l/2 . (14) 

Since a majority of extractions are unresolved point sources with very circular error-ellipses, 
we can ignore any correlations between uncertainties along each axis. To a good approx- 
imation, the uncertainty in radial separation between any two centroids can therefore be 
written: 

°r ^ \J°l+°l 3 ~ V / <+<> (15) 

where a^.) and (cr^ , a^.) are centroid variances in each axis for sources i and j re- 

spectively. A comparison between uncertainties in matched source radial separation (Equa- 
tion 15) and actual separations (Equation 14) after refinement is shown in the top panel of 
Figure 7 (for "10/20-match" case). They are both mutually consistent, although the spread 
is greater at separations < 0.3". After refinement and in the absence of systematics, any 
residual separation in a matched source pair must be due to intrinsic centroiding error alone. 



16 



Since separations between matched sources along each axis are to a good approximation in- 
dependently random and normally distributed with zero mean ((xi — Xj) ~ (Vi — yj) — 0), the 
quantity r (Equation 14) can be shown to follow a x distribution with two degrees of freedom 
(e.g., Evans et al. (2000)). This special case is also known as the Rayleigh distribution: 



where (3 is a parameter characterizing the width. This can be written in terms of the second 
moment (variance) of P(r) as follows 



For a given uncertainty ay as computed from Equation 15, the variation in the density 
of points with r along any horizontal cut in Figure 7 is qualitatively consistent with that 
predicted by Equations 16 and 17. 

The bottom panel of Figure 7 shows the dependence of the reduced x 2 (effectively the 
cost function in Equation 7 divided by number of degrees of freedom defined by Equation 8.) 
as a function of increasing number of images JVj in our 1000-image simulation. It's important 
to note that image offsets are not re-computed using repeated global minimizations for each 
new set of images JVj. Instead, the original full 1000-image solution of image offsets is used 
throughout to re-compute x 2 from Equation 7 as JVj is increased for all image pairs (to, n) 
such that n < to < JVj. As one approaches the full image set of JVj = 1000, one expects the 
reduced x 2 to converge to unity if, on average, residuals in source separations after refinement 
are purely consistent with extraction centroid uncertainties. The lower reduced x 2 values for 
smaller image numbers (and particularly for all image numbers in the "2/3- match" case) is 
due to the nonlinear behavior in x 2 when the number of degrees of freedom is small. Better 
fits (smaller x 2 ) can be obtained for very low numbers of degrees of freedom. In Figure 7 
(bottom panel), only the "10/20-match" case (with absolute refinement as studied above) 
shows approximate covergence to one, and the other curves are not far from it. Also shown 
is a case where only relative frame-to-frame matches and no astrometric references are used. 
This case tends to show a slightly higher reduced x 2 (— 1-04), which is significant since it is 
almost 20 standard deviations from the expected value in reduced x 2 (°x 2 = \/2/#dof, wnere 
#dof= 308935). This was traced as being due to slightly under-estimated undertainties in 
extraction centroids. This is not seen in the absolute refinement case (solid curve) since 
absolute astrometric uncertainties are themselves over-estimated and their (almost equal) 
contribution tends to lower the effective x 2 when combined with relative frame-to-frame 
matches. 




(16) 
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To summarize, we have presented a simulation to ascertain the performance of the 
pointingrefine algorithm. The model-dependent parameters entering our simulation can be 
isolated to properties of 3.6/im source populations, specifics of the IRAC band-1 array, 
such as PSF and pixelization, and a-priori telescope pointing knowledge. These can be 
appropriately rescaled to model other wavelengths and detectors. However, in the absence 
of systematics in the locations of potential matches between frames (both absolute and 
relative), and regardless of instrumental setup or detecter properties, the accuracy in refined 
pointing is purely dictated by the accuracy of point source centroids and match statistics. 
Our simulation (Figures 5 and 6) indeed shows that the refined pointing will typically never 
be worse off than the (combined) centroid uncertainties of matched pairs of sources. A well 
sampled and characterized PSF is expected to give centroiding accuracies to better than 
one-tenth of a resolution element. If errors are independently random, good match statistics 
can then only work in our favor to give the desired ~ 1/ V N matches improvement in pointing 
accuracy. 

4. Testing on Spitzer-IRAC Data 

The Infrared Array Camera (IRAC) is one of three focal plane instruments on the Spitzer 
Space Telescope (Fazio 2004). IRAC provides simultaneous ~ 5.2' x 5.2' images at 3.6, 4.5, 
5.8 and 8/im (bands 1-4). All four detector arrays in the camera are 256 x 256 pixels in 
size, with a pixel size ~ 1.2". We present here the results of a case study of observations 
acquired with IRAC during the in-orbit checkout period (Oct. 2003). In this section, we 
validate the pointing performance of IRAC and estimate the accuracy of refinement that can 
be achieved using a standard astrometric catalog and comparisons with our simulation of 
the 3.6/im band from Section 3. 

The observational request used for our case study consists of 105 regularly spaced images 
arranged in a rectangular raster pattern. The coverage map and geometry is shown in 
Figure 8. Adjacent images have ~ 20% overlap in each axis, with a coverage of 6 and 12 
pixels at the edges and (inner) corners respectively. The mapping was performed in repetitive 
horizontal scans as shown by arrows in Figure 8. The first image is at top right and the last 
at bottom left. All images across all bands were preprocessed for removal of instrumental 
signatures using the SSC's IRAC pipeline 11 and raw pointing information was attached to 
FITS headers. Source extraction was then performed using the SSC source extractor with 
PSFs characterized from in-flight data. Sources were extracted to a uniform signal-to-noise 



11 http: / / ssc.spitzer.caltech.edu/documents / SOM / 
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ratio of 5<r in each band, resulting in an average of ~ 39, 24, 13 and 6 extractions per image 
for bands 1, 2, 3 and 4 respectively. Errors in flux-weighted centroids were on average ~ 0.18, 
0.22, 0.26 and 0.27" (l-<x radial) for each band respectively. 

We used data from the Two-Micron All Sky Survey 12 (2MASS) to define a standard 
astrometric catalog for all IRAC bands. Zacharias et al. (2003) used the USNO CCD As- 
trograph Catalog (UCAC) (which is accurate to ~ 20 mas) to perform an assessment of the 
accuracy of 2MASS astrometry. They found random errors of 2MASS positions of ~ 85 to 
140 mas radial to a limiting K s magnitude of 15. We use astrometric sources detected in the 
2MASS K s band (2.12/xm) which is well suited to IRAC, and as a good working measure, 
we impose a magnitude limit of K s — 15. To this limit and at the highest galactic latitudes 

> 60°), one expects to find at least two astrometric sources per ~ 5.2' x 5.2' IRAC field 
~ 60% of the time. In the galactic plane anti-center, this increases to at least 30 sources 
80% of the time (J. Surace 2000, SSC, private communication). 

Using source extractions and astrometric references, the pointingrefine software was ex- 
ecuted on each of the four, band-dependent 105-image ensembles. Guided by a post-facto 
analysis of typical source separations and flux differences, point-source matching was per- 
formed using a nominal search radius of 2" and simultaneous flux matching performed with 
relative (frame-to-frame) and astrometric (frame-to-absolute) flux difference thresholds of 
4% and 50% respectively across all bands. Flux matching thresholds were set conserva- 
tively high to account for relative and absolute photometric calibration errors and instrinsic 
scatter between source populations detected in the 2MASS and IRAC bands. Relative and 
astrometric match statistics are summarized in Table 1. 

Figure 9 shows distributions in matched source radial separations before and after re- 
finement for all bands. Relative and astrometric matches have been separated. These distri- 
butions provide a powerful diagnostic with which to assess the in-flight pointing performance 
in the IRAC science instrument frame. The end-to-end pointing accuracy is a function of the 
inherent star-tracker accuracy, the spacecraft control system, how well the star-tracker bore- 
sight is known in the telescope pointing frame and focal plane array (science instrument) 
frame, and variations in these due to thermo- mechanical deflections. Spitzer's star-tracker 
assembly alone provides pointing and control to better than 0.3" absolute accuracy. The 
In-Orbit-Checkout (IOC) dataset analyzed here may have additional error sources because 
the images were obtained while the Spitzer pointing and data-acquisition systems were still 
being tuned and because the image integration times were very short (1.2 seconds). For ex- 
ample, if the estimate of the time needed for the spacecraft pointing system to "settle" were 
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slightly too short (i.e., underestimated), then the images would have been taken while the 
spacecraft pointing was changing slightly, providing an additional error term in the pointing 
transfer. More standard images with relatively long integration times would be much less 
affected by this source (as long as the settle time is small relative to the integration time). 
Therefore, the pointing accuracy we derive for this dataset may not be typical for Spitzer, 
but it still provides a valid test of the POINTINGREFINE algorithm. 

We can get an estimate of both the relative and absolute raw pointing for this IRAC 
dataset from the "before refinement" (dotted and thick lined) histograms in Figure 9. The 
main difference between bands is in the number of point source matches. Other than that, 
ranges in radial distributions are more or less consistent. Across all bands, the 1-a radial 
separation between frame-to-absolute matches is ~ 0.85 to 0.93" , while for relative (frame- 
to-frame) matches, this is ~ 0.61 to 0.65" (See Table 1). In addition to actual pointing 
dispersion, these estimates include a dispersion component from point source centroid errors. 
The contribution from centroiding error to the "before refinement" distributions however is 
negligible. Taking for instance the band-1 average centroiding error of ~ 0.18" and ~ 0.085" 
for 2MASS (1-a radial), this translates to effective match separation uncertainties of ~ 
0.14 and ~ 0.18" for frame-to-absolute and frame-to-frame (relative) matches respectively, 
assuming errors in each axis are uncorrelated (see Equation 15). This implies that the actual 
absolute and relative pointing for the frames in this observation is typically V0.85 2 — 0.14 2 ps 
0.83" and V0.61 2 — 0.18 2 « 0.58" (1-a radial) respectively. Similar results are found using 
other bands. We should emphasize that these estimates are strictly valid for this dataset 
alone. They are not representative of the actual pointing performance of IRAC where the 
absolute pointing is typically ~ 0". 45-0". 66 (1-a radial). For an in-depth study, see Barmby 
and Lacy (2004). 

The source separation distributions after refinement in Figure 9 allow us to validate how 
well (or whether) images were refined to within accuracies determined by source extraction 
centroids. As found in our simulation (Figure 5), images were refined to better than ~ 
65 mas (1-a radial) within truth pointings, with a corresponding ~ 200 mas dispersion 
in source separations after refinement (Figure 6). The top panel in Figure 10 shows the 
distribution of uncertainties in radial match separation (from Equation 15) as a function of 
separation r for band-1 after refinement. These are broadly consistent above the minimum 
cut-off uncertainty of ~ 100 mas imposed by the finite size of the PSF. The distribution 
in r at any given uncertainty cut a r can also be described by a Rayleigh distribution (see 
Equations 16 and 17). Furthermore, Figure 10 shows that typical systematic uncertainties 
due to inaccurately calibrated distortion are minimal, and if present, are expected to be 
much less than the centroiding errors. Guided by the simulation in Section 3 and Figure 9, 
we conservatively conclude that the majority of image pointings must be refined to better 
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~ 200 mas. In fact, we can predict the absolute dispersion in image pointings about truth 
using Equation 13. Given the typical match statistics listed in Table 1 and centroiding errors 
of 0.18" and 0.27" for bands 1 and 4 respectively, we expect dispersions of ~ 35 and ~ 140 
mas (1-cr radial) about truth for these bands respectively. We expect these dispersions to be 
smaller once centroiding errors are brought down using better characterized PSFs. 

In this case study, the number of astrometric (frame-to-absolute) matches are factors 
1.5 to 2 times greater than relative (frame-to-frame) matches (see Table 1). This is due to 
the specific flux difference thresholds used in source matching (4% and 50% for relative and 
absolute matches respectively; see above). The results of Figure 9 assumed 2MASS point 
sources with K s < 15. To ascertain the degree to which the number of absolute matches 
control the level of refinement, we repeated the source matching and refinement for band-1 
using a 2MASS magnitude cut of K s = 10. The bottom panel of Figure 10 shows source 
separations before and after refinement. The number of absolutes per frame was randomly 
distributed between and 4 (compared to ~ 20 relative matches/frame). When frame- 
to-absolute separations after refinement are considered alone (open histogram at bottom 
left), the mean separation is smaller by a factor of ~ 2 compared to the K s < 15 case 
(top left panel of Figure 9). This is due to inherently lower positional uncertainties for the 
brighter (K s < 10) sources. When combined with relative matches however, the mean source 
separation and hence level of refinement is essentially unchanged compared to our results 
above. 



5. Discussion and Conclusions 

We have presented a generic algorithm to perform astronomical image registration and 
pointing refinement. It is generic in the sense that it can be used on any set of astronomi- 
cal images which recognize the FITS and WCS pointing standards (Greisen and Calabretta 
2002; Calabretta and Greisen 2002). Either relative (self-consistent frame-to-frame registra- 
tion), single image absolute-astrometric, or simultaneous (relative and absolute) refinement 
is supported. The crux of the method involves matching point source positions between 
overlapping image frames and using this information to compute image offset corrections by 
globally minimizing a weighted sum of matched point-source positional differences. 

To ensure robust registration and refinement, the algorithm is best optimized with the 
following criteria. 

1. The random uncertainty in measured twist angle of an individual image frame (89) is 
assumed to be small such as to ensure sin 89 ~ 89. 89 < 60' is a good working measure 
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for the intended applications of this algorithm (where 1 — < 10 4 ). 

2. Input images have been accurately calibrated for distortion and possible nonuniform 
pixel scale. Any position-dependent systematic offset between source matches will limit 
the refinement accuracy to the size of the systematic error involved. 

3. Sufficient area overlap between adjacent image frames is needed to ensure good match 
statistics. 

4. Availablity of point sources with well defined flux distribution profiles approaching 
that of the instrument/detector's PSF. Extended sources will lead to larger centroiding 
errors. 

5. Well characterized PSF(s) for the image(s) at hand. These are crucial for accurate 
determination of source extraction centroids. If the inherent telescope pointing uncer- 
tainty is of order a third or larger than the detector pixel size, centroiding accuracies to 
better than one-tenth of a pixel or resolution element are recommended. For pointing 
uncertainties much less than the pixel size, there is little to be gained in resolution by 
improving the registration. 

6. With the suggested centroiding accuracy from step 5, at least five relative (frame- 
to-frame) and five absolute source matches per frame will give pointings refined to 
better than three- hundreths of a pixel (rms). The greater the number of matches, 
the better the refinement. A minimum of two source matches per image frame, either 
relative, absolute, or both is required to determine all offset parameters per image 
unambiguously. 

7. Astrometric (absolute) reference sources should be used wherever possible. These 
provide a baseline to counteract any systematic deviations from "truth" , or expected 
pointing in the ICRS, especially if absolute refinement is desired. 

8. It is assumed that uncertainties in image pointing and point source extraction centroids 
are random and independent. Prior image pointing uncertainties should be used if 
available. These will prevent from erroneously over-correcting the pointing for cases in 
which it is well known a-priori. 

Our simulations show that potentially good refinement can be obtained with minimal 
requirements. For a large fraction of images in an ensemble, refinements of better than 
~ 65 mas (1-a radial from truth) can be obtained with an average of ~ 10 relative and 
~ 20 absolute matches per frame with extraction centroids < 0.2". This amounts to an 
improvement of 95% relative to truth for a majority of images. This could be better with 
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higher match statistics and/or more accurate centroids, since typically the 1-a dispersion 
about truth scales as ~ a ext ^2/N ext for a given centroiding uncertainty a ext and number of 
matches N ext . 

Analysis of observations from Spitzer's IRAC instrument shows that the dispersion in 
source separations after refinement is entirely consistent with the inherent dispersion in 
extraction centroid uncertainties. This implies that systematic uncertainties such as inac- 
curately calibrated distortions are negligible, since otherwise dispersions in matched source 
separations after refinement would be larger relative to centroiding errors. Comparing dis- 
persions of refined pointings about truth with matched source separations as found in our 
simulation, and rescaling to the appropriate numbers of matches using 1 / V N match scaling, we 
predict (at the 2a level) refinements to better than ~ 70 mas and ~ 280 mas for IRAC bands 
1 and 4 respectively. These bands bracket two extremes in available source matches, and 
these refinement estimates correspond to ~ 55 and ~ 8 (astrometric and relative) matches 
per band-dependent frame respectively. 

The goal of astronomical image registration is to exploit the resolution capabilities of 
existing and upcoming detectors whose pointing control and stability may not evolve at 
the same rate. The aim is to optimize the achievable signal-to-noise and science return 
therein. The algorithm presented herein is just the tip of the iceberg at exploring one of 
many optimization techniques used in the diverse fields of image and signal processing and 
computer vision science. 
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makes use of data products from the Two Micron All Sky Survey, which is a joint project of 
the University of Massachusetts and the Infrared Processing and Analysis Center/California 
Institute of Technology, funded by the National Aeronautics and Space Administration and 
the National Science Foundation. 



A. Elements of the Coefficient Matrix 

In this section, we provide general analytic expressions for elements of the coefficient 
matrix M (Equation 10). These elements are obtained by applying the minimization condi- 
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tions (Equation 9) to the cost function denned by Equation 7. The "base" coefficent labels 
A, B and C correspond to the three equations obtained by evaluating the partial derivatives 
in Equation 9 (labelled equations A, B and C from left to right respectively). Position and 
uncertainty variables appearing in the expressions below were defined in Section 2.2. 
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The Error-Covariance Matrix 



The full error- covariance matrix is one of the by-products of the pointingrefine software. 
This reports all variances and covariances for and between all (inter and intra) image offsets 
necessary for refinement from the global minimization. It can be used to explore the strength 
of long-distance correlations between images in a mosaic and the presence of any undue 
systematic walks after refinement. The latter could arise if one lacks the desired number of 
absolute astrometric sources, or correct magnitude of prior image pointing uncertainties as 
was discussed in Section 2.3. 

For a given pair of images (i, j) and three computed offsets per image (59, 5X, 5Y), 
we have a total of six possible (correlated) offset pairs or covariance matrices. The possible 
covariances (or variances if % = j for the same offsets) for any two images are therefore 



cov(9 i ,9 j ), cov(9 i ,X j ), cov(9i,Yj), cov(X h Xj), cov(X h Yj), cov(Y h Yj 



(Bl) 



If we define any of these covariances generically as cov(ai, f3j) then the format for the error- 
covariance matrix for all possible image pair combinations (i, j) where (i = l,2...m; j = 
1, 2...m) is as follows: 



/ cov(ai,Pi) cov(a 1 ,/3 2 ) 
cov(a 2 ,/3i) cov(a 2 ,/3 2 ) 



cov(ai, f3j) 



. cov(ai,(3 m ) \ 
. cov(a 2 ,p m ) 



\ cov(a m ,(3i) cov(a m ,p 2 ) 



(B2) 



cov(a m ,p m ) J 



Figure 11 shows greyscale representations of covariance matrices for all combinations of 
offset parameters (Equation Bl) for the 105-image IRAC test case assuming both relative and 
absolute astrometric source matches in the refinement. With 105 images, each covariance 
matrix has 105 x 105 elements (or pixels in this representation). It is important to note 
that the cross-correlation between different offset types is not symmetric. For instance, the 
correlation between 9-i and Xj is not the same as 9j and Aj. Figure 12 shows the same set but 
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with only relative (frame-to- frame) matches used for the refinement. Comparing Figures 11 
and 12, we draw the following conclusions: first, for the case including astrometric matches, 
"long-distance" correlations in image offsets are greatly reduced. This is due to astrometric 
sources anchoring each image to the fiducial reference frame, making them more indepen- 
dent of each other. The covariance matrices become essentially block-diagonal. Second, 
uncertainties (variances along the diagonal for the same offset pair combination) are greatly 
attenuated when astrometric matches are used compared to the relative-only match case. 
There are greater numbers of degrees of freedom per image when astrometrics are included, 
and this reduces the relative uncertainty. 

The relative-only case exhibits greater long-distance correlations (larger off-diagonal 
values) since image positions are dictated solely by frame-to-frame matches in image over- 
lap regions. This makes each successive image position dependent on it's nearest neighbor 
positions, which depend on their own neighbours and so on throughout the system of linked 
images. Apart from astrometric matches reducing long distance correlations, this can also 
happen if prior image pointing uncertainties are intrinsically smaller. Small priors will pull 
refinement offsets from the global minimization towards zero, which minimizes the L apr i or i 
term in Equation 4. The priors force a constraint on each individual image to prevent "over- 
refinement", regardless of the number of source matches present. Each image therefore 
becomes more independent of its neighbors, and long-distance correlations are reduced. 

We also note the rich and diverse patterns in the covariance matrices for this observation, 
especially the relative-only case (Figure 12). These are characteristic of the image layout and 
mosaic map geometry (see Figure 8). For instance, take the cov(Yi,Yj) matrix in Figure 12. 
The checker-board pattern arises from the relative image numbering in the map and how 
this translates to the numbering of elements in the covariance matrix (Equation B2). In 
the mosaic, the images repeat from right to left to create a leg, then left to right and down 
again, for a total of seven legs. Adjacent image pairs (i, j) along vertical sections in the 
map are strongly correlated in their (Y^, Yj) offsets, while widely separated images are less 
correlated. For instance, image number one at top right has its Y offset strongly correlated 
with the Y offsets of images directly below it (i.e., image numbers 30, 31, 60, 61, 90 and 
91). The high correlations are therefore with every ~30th image giving bright regions in the 
greyscale covariance image. Also, the Y offset of image number one is least correlated with 
that of images in the far left vertical strip (image numbers 15, 16, 45, 46, 75, 76 and 105), 
resulting in dark regions in the covariance image. 

In the end, one purpose of the covariance matrix is to visualize the degree of correlation 
between image positions in the ensemble as a whole. This allows one to ascertain whether 
refinement solutions are driven by any particular invalid input assumptions (e.g., priors), or 
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insufficient astrometric reference source information if robust absolute refinement is desired. 
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Inputs: FITS image data, point source 
extraction tables, optional table of 
astrometric source references and 
software control parameters. 



If and only if input source positions 
are given in terms of "undistorted" x,y 
pixel positions, convert these to RA, 
Dec to perform source matching. 



Perform point-source position and 
flux matching between all possible 
overlapping pairs of image frames in 
the RA, Dec system (including 
astrometric references if specified). 



Define reference image frame for 
registration and refinement: 
If absolute references specified, use 
bounding box defined by all image 
inputs containing astrometrics, 
otherwise, use specific single input 
image, or, maximally correlated image 
found by software to perform relative 
refinement. 



Transform correlated point source 
positions and uncertainties to tangent 
projection plane defined by reference- 
image frame. 



Transform a-priori image pointing 
uncertainties/twists to equivalent 
reference image frame measures. 



Compute matrix elements for global 
minimization computation and setup 
linear sparse matrix solver. 



Compute translational and rotational 
offsets of all image inputs in reference 
image frame using matrix solver. 



If matrix is singular, go back 
and use only frame-to-absolute 
matches on a single frame 
basis if absolutes specified, 
otherwise abort. 



Compute error-covariance matrix for 
translational and rotational offsets 
using matrix solver. 



Use translational/rotational Cartesian 
offsets to correct for input image 
positions/orientations in reference 
image frame. 



Transform corrected Cartesian 
positions/orientations of images to 
equivalent refined celestial pointings, 
position angles and CD matrix 
elements in the WCS standard. 







Convert error-covariance matrix for 
Cartesian offset corrections to 
uncertainties in refined WCS. 







Outputs: write refined pointing 
keywords and CD matrix elements to 
FITS headers of input images, 
optionally output tables of refined 
WCS parameters, Cartesian offset 
corrections, point source matches, 
refinement statistics and diagnostics. 



Fig. 1. — Processing and algorithmic flow in pointingrefine software. 
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Fig. 2. — A simple three-image mosaic. Filled circles are sources detected in image n and 
open circles are sources in m or 1. 
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Fig. 3. — One-dimensional representation of projection geometry showing images on sky and 
in tangent plane of reference image. Images 1 and 2 have the same physical size but different 
projected sizes. 
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Fig. 4. — Unrefined mosaic section of 1000-image IRAC simulation at 3.6yum on left and 
same section after pointing refinement on right. Field dimensions are ~ 2.8' x 3.8'. 
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Fig. 5. — Top: distributions of image center separations relative to "truth" before and after 
refinement. Bottom: Offsets in RA and Dec between image centers relative to truth. Circle 
represents 2-a region of truth - refined distribution for case with ~ 2 relative/3 absolute 
matches per image (dot-filled histogram in top figure). 
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Fig. 6. — Top: distributions of the magnitude of refinement, represented as a percentage of 
initial truth - unrefined image separation. Bottom: distributions of matched source radial 
separations before and after refinement. 
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Fig. 7. — Top: uncertainty in matched source radial separation (from centroid uncertain- 
ties) as a function of actual source separation after refinement. Bottom: reduced \ 2 (i- e -> 
X 2 /number of degrees of freedom) as a function of mosaic subset composed of successively 
increasing numbers of images % (see Section 3). 
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Fig. 8. — Coverage mosaic and mapping geometry of the 105-image IRAC observation used 
in this analysis. Adjacent images have ~ 20% overlap with a coverage of 6 and 12 pixels at 
edges and (inner) corners respectively. Mapping direction is shown by arrows starting at top 
right and ending at bottom left. Open circles are image centers. 
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Fig. 9. — Distributions of matched source radial separations before and after refinement for 
all bands of Spitzer's IRAC instrument. Relative and absolute astrometric (with 2MASS 
magnitudes K s < 15) matches have been separated. 
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Fig. 10. — Top: uncertainty in matched source radial separation (from centroid uncertainties 
of matches) as a function of actual source separation after refinement for IRAC band-1. 
Bottom: same as Figure 9 for IRAC band-1 but with astrometric matches having a 2MASS 
magnitude limit of K s = 10. 




Fig. 11. — Greyscale representation of covariance matrices for all orthogonal and rotational 
offsets in IRAC test case when absolute astrometric matches are used. See Equation B2 for 
matrix definition. 
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Fig. 12. — Greyscale representation of covariance matrices for all orthogonal and rotational 
offsets in IRAC test case when absolute astrometric matches are not used (i.e. only relative 
frame-to- frame matches used). See Equation B2 for matrix definition. 
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Fig. 1. — Processing and algorithmic flow in pointingrefine software. 

Fig. 2. — A simple three-image mosaic. Filled circles are sources detected from image n and 
open circles are sources from m or 1. 

Fig. 3. — One-dimensional representation of projection geometry showing images on sky and 
in tangent plane of reference image. Images 1 and 2 have the same physical size but different 
projectedsizes. 

Fig. 4. — Unrefined mosaic section of 1000-image IRAC simulation at 3.6yum on left and 
same section after pointing refinement on right. Field dimensions are ~ 2.8' x 3.8'. 

Fig. 5. — Top: distributions of image center separations relative to "truth" before and after 
refinement. Bottom: Offsets in RA and Dec between image centers relative to truth. Circle 
represents 2-a region of truth - refined distribution for case with ~ 2 relative/3 absolute 
matches per image (dot-filled histogram in top figure). 

Fig. 6. — Top: distributions of the magnitude of refinement, represented as a percentage of 
initial truth - unrefined image separation. Bottom: distributions of matched source radial 
separations before and after refinement. 

Fig. 7. — Top: uncertainty in matched source radial separation (from centroid uncertain- 
ties) as a function of actual source separation after refinement. Bottom: reduced x 2 (i- e -> 
X 2 /number of degrees of freedom) as a function of mosaic subset composed of successively 
increasing numbers of images % (see Section 3). 
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Fig. 8. — Coverage mosaic and mapping geometry of the 105-image IRAC observation used 
in this analysis. Adjacent images have ~ 20% overlap with a coverage of 6 and 12 pixels at 
edges and (inner) corners respectively. Mapping direction is shown by arrows starting at top 
right and ending at bottom left. Open circles are image centers. 

Fig. 9. — Distributions of matched source radial separations before and after refinement for 
all bands of Spitzer's IRAC instrument. Relative and absolute astrometric (with 2MASS 
magnitudes K s < 15) matches have been separated. 

Fig. 10. — Top: uncertainty in matched source radial separation (from centroid uncertainties 
of matches) as a function of actual source separation after refinement for IRAC band-1. 
Bottom: same as Figure 9 for IRAC band-1 but with astrometric matches having a 2MASS 
magnitude limit of K s = 10. 

Fig. 11. — Greyscale representation of covariance matrices for all orthogonal and rotational 
offsets in IRAC test case when absolute astrometric matches are used. See Equation B2 for 
matrix definition. 

Fig. 12. — Greyscale representation of covariance matrices for all orthogonal and rotational 
offsets in IRAC test case when absolute astrometric matches are not used (i.e. only relative 
frame-to- frame matches used). See Equation B2 for matrix definition. 
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Table 1. Statistics for IRAC-observation case study. 



A-Band 


(#Afo.) a 




(D bef ) A c 


(D bef ) R c 


(D aft ) T d 


x 2 


doP 


X 2 /dof 


fim 






arcsec 


arcsec 


arcsec 








3.6 


34.7 


20.5 


0.847 


0.605 


0.151 


8580.84 


10821 


0.792 


4.5 


20.7 


13.7 


0.883 


0.622 


0.152 


3843.93 


5645 


0.681 


5.8 


9.4 


7.6 


0.933 


0.603 


0.142 


1250.31 


2443 


0.511 


8.0 


5.3 


2.1 


0.913 


0.651 


0.203 


1240.50 


1249 


0.993 



a Average number of absolute matches per frame. 

b Average number of relative (frame-to-frame) matches per frame. 

c (Dbef)A and (Dbef)R represent the mean matched source separation before refinement for 
absolute and relative matches respectively (see Figure 9). 

d (D a f t ) T represents the mean total matched source separation (for both relative and abolute 
matches) after refinement. 

e Number of degrees of freedom (see Equation 8). 



